New Multi- order exact solutions for a 
class of nonlinear evolution equations 



Bijan Bagchi^* Supratim Das^^and Asish Ganguly^^ 

''^Department of Applied Mathematics, University of Calcutta, 
CN : 92 Acharya Prafulla Chandra Road, Kolkata-700009, India 

■ Department of Mathematics, Indian Institute of Technology, Kharagpur 



(N 



Kharagpur-721302, India 



Keywords: Travelling waves; Lame equation; Lame function; Jacobi elliptic function; 
Elliptic integral function. 



C/0 : Abstract 

We seek multi-order exact solutions of a generalized shallow water wave 
equation along with those corresponding to a class of nonlinear systems de- 
scribed by the KdV, modified KdV, Boussinesq, Klein-Gordon and modified 
^SJ ■ Benjamin-Bona-Mahony equation. We employ a modified version of a general- 

^ . ized Lame equation and subject it to a perturbative treatment identifying the 

solutions order by order in terms of Jacobi elliptic functions. Our multi-order 
\^ • exact solutions are new and hold the key feature that they are expressible in 

. terms of an auxiliary function / in a generic way. For appropriate choices of / 

we recover the previous results reported in the literature. 
PACS number (s): 02.30. Jr, 02.30.Ik, 05.45.-a 

^ 1 Introduction 

d \ Seeking tractable solutions of nonlinear evolution equations has been the focus of 

intense study for the past several decades (see, e.g. , [lj-[II]). While many strategies 
have been employed, such as the Hirota bilinear method [2], homogeneous balance 
method [3], trigonometric method [1], the hyperbolic method [5], the Jacobi elliptic 
function method [B] and iterative method [HIS] to name a few, perturbative techniques 
have also been used to extract multi-order exact periodic solutions based on Lame 
equation and Jacobi elliptic functions [9], [TOl [TT] . In this note we propose a generalized 
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class of Lame equation comprising an auxiliary function f{^) to consider correlations 
between it and the perurbatively reduced nonlinear evolution equation for different 
orders of the perturbative parameter. 

A generalized Lame equation is expressed as 

fl+[X-(n + l){n + 2)aMfiOr] V = 0, (1) 

where A stands for an eigenvalue, n is a positive integer and we restrict the auxiliary 
function /(^) to obey an elliptic equation 

(|)'=(l + «i/^)M^ + «3). (2) 

In ([T]) and ([2]), ^i, 0,2, 03 are real constants. It is well known that equation ([2]) admits 
of several categories of solutions for different values of ai, 02, as all expressible in 
terms of modulus k of the Jacobi elliptic functions. For later use we summarize a few 
relevant ones in Table [T] 

Table 1: /(^) is provided along with the values of the parameters aj ,j = l, 2, 3. The 
complementary modulus is denoted by k'"^ = 1 — k"^ |12j . 
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sn(^, k) 

cn{i, k) 

dn{^, k) 
nc(e, A;) I.e. 
ndiC, k) I.e. ^^1^^^^ 
sc(e, k) I.e. J^^J 
cd{i, k) I.e. 
dc{i, k) I.e. ^^1'^ 
sd{i, k) I.e. ^ jl^j 



It should be mentioned that while we would attend to the enlisted set of /(^) 
in Table [1] to widen the scope of enquiry, some particular cases of f{C)i namely 
sn{^, k), dn{^, k) and cd{^, k), have been studied in the literature [9l[T0l[Tl] to arrive 
at solutions of the relevant PDEs by applying the perturbation method on the latter. 
However, as is evident from Table [T] there do remain other variants of /(^) which 
open up different cases of new solutions not only for the already considered PDEs 
but other types as well by perturbatively reducing them to an ODE form so as to 
match with ([1]). 

We have selected a physically meaningful model of a generalized shallow water 
wave equation [131 [HI [13 [IS] which is amenable to a perturbative treatment by look- 
ing for a travelling wave variable u{x, t) = u{^) with ^ = j{x — ct) where 7 is the wave 
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number and c is the wave speed. We discuss the procedure of obtaining the solutions 
in the next section. In the following section we also take up applicability of our 
scheme to other nonlinear evolution equations namely, the Korteweg de Vries (KdV) 
equation, modified KdV equation, Boussinesq equation, Klein-Gordon equation and 
modified Benjamin-Bona-Mahony (mBBM) equation with a view to tracking down 
new hitherto unexplored multi-order solutions. The main point of our analysis is to 
demonstrate that all such solutions can be written down in a generic form in terms 
of the function /. 

To make equation ([T]) accessible in terms of known Lame functions, we recast it 
in terms of a new variable t] by applying the transformation /(^) = y^. Using 
this gives 



cPy 1 



1 1 1 

- + + 

r] 11+- r] + fi 



dy _ u + {n+l){n + 2)r] ^ ^ 



where the parameters u and u are defined by u = — and u = — —. 

^ " " a2 aia2 

Equation ([3]) is readily solvable for the special cases of n = 1 and 2 for different 
choices of fi and u subject to certain relation between them. The solutions expressible 
by the corresponding Lame functions are [9] 



n 



L{{0 = (l + ai7/)^(a3 + a2?7)^ = (1 + aif )^(a3 + )^ (4) 



A = — (a2 -|- aitta) [provided u = -\ )]. 

ai 



n = 1 L{^(e)=?72(i + ai?7)2 =/(l + ai/^)2, (5) 
A = — (a2 + 4aia3) [provided z/ = (4/i H )]. 



n = 1 Li''{0 = vHa, + a2r]r^=f{a3 + a2fr^, (6) 

4 

A = — (4a2 + flitts) [provided = (/i H )]. 



n = 2 L2(0 =?7^(l + ai77)5(a3 + a277)^ =/(l + ai/2)^(a3 + a2/')i (7) 
A = — 4(a2 + cuo-s) [provided u = 4(yU H )]. 

Henceforth we will be guided by the solutions ([4]) - ([7]) along with an appropriate 
/ from Table 1 for the various multi-order cases that follow from a nonlinear PDE. 
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2 A generalized shallow water wave equation 



We turn attention to a generalized class of shallow water wave (GSWW) equation 
given by [ISl H [151 [IS] 

Uxxxt + aur^u^t + (^UtUr,^ - u^t - u^x = 0, (8) 

where a, /3 G M — {0}. Under Boussinesq approximation the derivation of ([H]) results 
from the classical study of water waves. More interestingly, there also follows various 
classical and non-classical reductions from GSWW such as the KdV and BBM equa- 
tions. Investigations of Painleve tests reveal complete integrability for specific values 
of a and /3 namely Q; = /3ora = 2/3 [13]. Recently we extended [TB] the Jacobian 
elliptic function method to classify new exact travelling wave solutions expressible in 
terms of quasi-periodic elliptic integral function and doubly-periodic Jacobian elliptic 
functions. 

For the travelling wave solutions u = u{^), equation (|8]) can be reduced in terms 
of the variable v{^) = ^ to the ODE form 



^ + Pv^ + Qv + ci = 0. (9) 

where P = Q = ^ and ci is an integrating constant. 
To tackle ([9]) perturbatively we set 

v{0 = vo{0 + evi{0 + e'v2{0 + --- (10) 

where e(> 0) is a small parameter and fo(0) "^1(0 5 ''^2{0i ■■■ represent various multi- 
order solutions like the zeroth-order, first-order, second-order solutions etc. of equa- 
tion Accordingly we can write u{^) as 

«(0 = %(e) + e«i(0 + eV(0 + --- (11) 

where Ui{^) = Vi{T) dr. 

Substituting the series (|T0|) in equation ([9]) we obtain for each power of e the 
corresponding equations 



^ • -77^ + P% + Qvo + ci = 0, (12) 



t2 

: -^ + {2Pvo + Q)vi = 0, (13) 



: ^ + {2Pvo + Q)v, = -Pvl (14) 



and so on. 
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We now proceed to solve the above chain of equations by first expanding vq namely, 

I 

V, = Y,Arf\ (15) 

i=0 

where Aj's are constants and then comparing highest order linear and nonlinear terms 
in f lT^ . In this way we obtain / = 2 by making use of (j^]). Thus f lT5]) gets reduced to 
the quadratic form 

v^ = A^ + A^f + A2f. (16) 



2.1 Zeroth-order exact solution 

With Vq given by ( IT6|) we are led to the following system of coupled equations: 

2a^A2 + PAl + QAQ + ci = 0, (17) 

A,[{a2 + aia3) + 2PAo + Q] = 0, (18) 

4:{a2 + aia3)A2 + P{Al + 2AoA2) + QA2 = 0, (19) 

A^[aia2 + PA2] = 0, (20) 

A2[6aia2 + PA2] = 0. (21) 

It is readily seen that we generate the following set of consistent solutions Aq = 
— ^[Q + 4(a2 + aifls)], Ai = and A2 = — ^"p"^ along with the constraint 

[16(a2 + aios)^ - 48010203] - + ^Pci = 0. (22) 



From fl22l) we obtain wave speed as 



[1 ± 472{(a2 + oiag)' - 3010203 + ^^c,p]-' (23) 

07 



can always adjust the integration constant ci to keep (02 + 0103)^ + ^^ci > 8010203. 
Further for the finiteness of c we require (02 + OiOs)^ + ^^Ci 7^ 8010203 + 
From ( 1121) and flTB]) we get for the first integral 

1 M 6O1O2 „2 /r^.N 

^0 = ^ = [Q + 4(02 + 0103)] - -^f (24) 
which in turn gives the zeroth-order solution 

no = -^[Q + 4(a2 + ai03)]e -^jf (25) 

As is evident, uo depends upon the choice of the auxiliary function /. In Table [2] 
we furnish the various forms for uq corresponding to different elliptic functions of 
Table [H A sample graph of the zeroth-order solution uq for / = sn{^, k) is depicted 
in Figure [T] and a plot of wave-speed for variation of the constant Ci is shown in 
Figure |21 
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Figure 1: Zeroth-order exact solution corresponding to f{^) = sn{^,k) for k = 
0.5, a = 1, 13 = A, 7=1 



Figure 2: Wave speed c vs constant of integration Ci for f{^) = sn{^, k) for k = 
0.5, a = 1, /3 = 4, 7 = 1 

2.2 First-order exact solution 

Knowing vq from fl2ll) . equation f|T3l) reads 

^ + [-4(a2 + aiaa) - 12a^a2f]vl = 0, (26) 

which matches with ([1]) for A = — 4(a2 + aifls). This gives n = 2 implying that the 
solution of fl2^ can be written as 

vi = ^ = = C2/(l + a.nkl + -f)K (27) 

where C2 is an arbitrary constant. 

Integration of (127|) with the use of ([2]) immediately provides 

ui = ±|/^ (28) 



Table 2: Wave velocity c = , , i?(</>, k) is the incomplete elliptic integral 

function of second kind where sinci) = sn{^, k) and A. = —-^[Q — 4(1 + k"^)]- 





Zeroth-order solution uq 


sn(^, A;) or cn(^, A;) or dn{^, k) 
nc{(,, k) or sc(^, A;) or dc{^, k) 
nd{^, k) or A;) or cd{^, k) 


A^ + ^E{(l),k) 
AC + j,[E{cl), k) - k)dc{i, k)] 
Ae + |[^((^, A;) - k^sn{i, k)cd{C, k)] 



which serves as a first order approximation to equation ([8]) . Figure |3] gives a plot of 
ui for / = sn(^, A:). 




Figure 3: First-order exact solution corresponding to f{^) = sn{^, k) for k = 0.5, a 

1, /3 = 4, 7 = 1, C2 = l 



2.3 Second-order exact solution 

Using (l24l) and (127|) the second-order equation (fT4|) takes the form 



[-4(a2 + aiag) - 12aia2/1t'2 = -Pc^/ll + aj'){as + a^f 



(29) 



We observe that the coefficients of (1291) are polynomials in /. This prompts us to 
take a polynomial ansatz for V2 namely 



V2 



(30) 



1=0 



Substitution of the above into we get / = 4 on equating the highest power of / 
from both sides. Thus (I5U]) is reduced to the biquadratic form 

V2 = B^ + B^f + B^f + Bs/' + B,f\ (31) 
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Putting ( 13T]1 into fl29ll and equating the coefficients of /* (z = to 6) to zero, we 
have the following system of coupled equations, 



a3i?2 - 2(a2 + aia3)5o 
(as + aia^)Bi - 203^3 
1203^4 — 12aia2-Bo + -Pc2a3 
(a2 + aia3)i?3 - 2aia2i?i 
-6aia2-B2 + 12(a2 + aia3)-B4 + -Pc2(a2 + aiQs) 

aia2i?3 
8aia2i?4 + Pc2aia2 

Solving for So, -Bi, -B2, -B3, -B4 we obtain 



0, 
0, 
0, 
0, 
0, 
0, 
0. 



a-jclP 
24aia2 



, Si = 0, S2 



[a2 + 0103)02^ 
12aia2 



, ^3 = 0, S4 



(32) 



(33) 




Figure 4: Second-order exact solution corresponding to f{^) = sn{^, k) for k 
0.5, a = 1, /3 = 4, 7 = 1, C2 = 1 



Hence the solution of equation fl29l) can be written as 



V2 



24aia2 

Using ([2]) we can integrate flM|) to get 



[03 + 2(a2 + aifls)/^ + 3aia2/^ 



U2 



24aia2 

^2 is plotted for / = sn{^, k) in Figure IH 



(34) 



(35) 



8 



3 Application to other nonlinear 
evolution equations 

In the previous section, we considered a nonhnear second-order ODE containing 
a functional parameter / and obtained specific solutions in terms of /. It was shown 
that the use of Jacobi elliptic function method in a perturbative way yielded new 
multi-order exact solutions that depended upon / in a generic way. In this section 
we show that our approach can be profitably applied to other nonlinear equations 
resulting in more general types of solutions that have not been reported to the best 
of our knowledge. A complete list of our results are presented in Table [31 



Table 3: New general multi-order exact solutions for various non-linear evolution 
equations 



Equation 


Multi-order exact solutions m„, n = 0, 1, 2 


(z) KdV equation 

Ut + auUx + (iuxxx = 


Wo = f - ^[(^2 + aias) + 3aia2p], 
[where c = ±(3p{16{a2 + aia-sY — 48aia2a3 + 2^ci}2] 
ui = C2/(l + ai/2)5(a3 + a2p)^ 

2 

C'XCX. V (~\ / \ f'? o 

^2 = 48aia2/372 + 2(a2 + aiag)/^ + 3aia2/*J 


(n) mKdV equation 

Ut + av?Ux + Puxxx = 


[where c = /3p{a2 + aia^)] 
Ml = C2(l + ai/2)3(a3 + a2p)^ 

«2 - ± 1 y 6/37"aia2 /[(«2 + Olfls) + 20102/^] 


(in) Boussinesq equation 

'^tt CqUxx Ol'^xxxx )xx ~ 


Uo = -^[cq - + 4072(02 + 0103)] - ^aia2p, 
[where c = {cq ± Aap^y{a2 + 0103)2 — 3010203} 2] 
ui = C2/(l + oi/2)i(a3 + 02/2)^ 

^2 = 24aia2a72 [03 + 2(o2 + ai03)/2 + 30ia2/^] 


(iv) Klein-Gordon equation 

Utt - Uxx + au + /3u^ = 


kherec=±{l ^,.(,,:„,„3) }^] 

Ml = C2(l + Oi/2) 2 (03 + 02/2) 2 
^2 - ± 1 \/ 2(c2-l)72aia2 /[(«2 + OlOs) + 20i02/2] 


(v) mBBM equation 

Ut + CqUx + U^Ux + fSUxxt = 


Uo = ±^6/37200102 /, 

[^he^eC= l+(a2+aia3w] 
Ml = C2(l + 01/2)5(03 + 02/2)5 

U2-±^^ 6/372Lia2 /[(«2 + ^103) + 20i02/2] 
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(z) KdV equation 



The concerned equation is of the form 

ut + auucc + I3uxxx = (36) 

where a and (3 are real parameters. Let us choose / = cd{S^,k) for which ai = —1, 02 = 
— fc^ and 03 = 1. Substituting these values into the zeroth-order, first-order and second- 
order exact solutions listed in Table 3, we find 

no = ^ + '-^{l + k')-'-^k^cd\(,k) 
a a a 

m = C2k'^cd{C,k)sd{C,k)nd{C,k) (37) 
= -^^^^[l-2(l + A;2)cd2(^,A;)+3A;2cd4(e,A:)]. 



These solutions have been obtained in [10] for a = 1 where the authors used the notations 
m and k in places of k and 7 respectively and their scale is A = k''^C2, ci = 0. The 
KdV equation has also been studied in [S] and their results can be extracted from ours for 
/ = sn(^, k) corresponding to a = 1. 



(ii) mKdV equation 



Ut + au'^Ux + /3uxxx = 0. (38) 



Our general result listed in Table [3] contains, as special cases, the following results reported 
in earlier works. 

• / = sn{C,k) (Ref [9j) 



Uo = ±-fk\l-— sn{^,k) 
V a 

Ul = C2 cn{S^, k)dn{^, k) (39) 
(l + k'^)cl I 6a 2A;2 ,n -,n 



/ = cd{i, k) (Ref [in 



6 /3 

Wo = ±7^a/ cd{£,,k) 

V a 

= {I - k^)c2 sd{i,k)nd{i,k) (40) 

U2 = iS^^V-y^'^trn^^'^^'"^-'] 
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{in) Boussinesq equation 

utt - clu^x - auxxxx - /3iu^)xx = 0. (41) 

From the general solution given in terms of / (see Table [3]) one can deduce the ones obtained 
in [S] as special case for / = sn(^, k). 



(iv) Klein-Gordon equation 

Utt - Uxx + au + fiu^ = 0. (42) 

The above equation was considered in [11]. Their results correspond to the case / = dn(^, k) 
and we do indeed recover their solutions from Table 3 as given below: 



-0 = W ^^2(p_2) ^"(^'^^ 



Ul = k^C2 sn{^, k)cn{^, k) (43) 



U2 = T§J^;^^T^^—^^dn{tk)[2-k'-2dn\tk)] 



(v) mBBM equation 



Ut + CqUx + V?Ux + I^Uxxt = 0. (44) 
For / = sn{^,k) the multi-order solutions to ()44p read from Table 3 

uq = ib7Y^ 6c/3 sn(^, k) 

Ul = C2 cn{S,, k)dn{S,, k) (45) 
These results were obtained in Ref [9]. 



4 Concluding remarks 

We have reported exact multi-order periodic solutions for a generalized shallow wa- 
ter wave equation based on various types of Jacobi elliptic functions in a perturbative 
framework. The presence of an auxiliary function, which makes reference to any particular 
Jacobi function explicit, can be well adjusted to connect with other classes of PDE such as 
the KdV, modified KdV, Boussinesq, Klein-Gordon and modified Benjamin-Bona-Mahony 
equation. The general character of our results encompasses those special cases which have 
earlier been studied in the literature. 
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